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Abstract. - It is shown how S-matrix theory and the concept of continuous quantum measure- 
ments can be combined to yield Markovian master equations which describe the environmental 
interaction non-perturbatively. The method is then applied to obtain the master equation for 
the effects of a gas on the internal dynamics of an immobile complex quantum system, such as a 
trapped molecule, in terms of the exact multi-channel scattering amplitudes. 
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C| Introduction. — A truism of quantum physics tells 
Q4hat no system is perfectly isolated and it is therefore 
+— > not surprising that the study of open quantum systems 
^ is an ubiquitous theme of present day quantum mechan- 
^ics, see [1-4] and refs. therein. An important class of 
Q-fevolution equations for open systems are Markovian mas- 
L^ter equations. They imply that environmental correla- 
. tions disperse fast, so that on a coarse-grained timescale 
^ the temporal change of the system state p depends on the 
^ present state of the system, but not on its history. From 
the strict point of view of an operationalist (who dismisses 
the notion of a "system" altogether and takes p as describ- 
ing an equivalence class of preparations in the lab [5]) one 
may even argue that any valid differential equation for the 
time evolution of p must generate a completely positive 
dynamical semigroup [6] and must hence be Markovian. 

Putting the pros and cons of Markovian vs. non- 
Markovian formulations aside, it is fair to say that a large 
class of open quantum systems is described appropriately 
by time-local master equations. At the same time, it is cu- 
rious that the Markov property does not emerge naturally 
in standard microscopic derivations. Rather, one has to 
impose it "by hand" , usually by interpreting some quanti- 
ties as correlation functions, which must then be assumed 
to be (5-correlatcd. This may be still transparent in weak 
coupling calculations such as the Bloch-Redfield approach 
[1], but tends to be awkward if a non-perturbative treat- 
ment of the interaction with the environment is needed. 

In the present letter I would like to motivate and ex- 
emplify a general method of obtaining master equations 



which do incorporate the microscopic interactions in a 
non-perturbative fashion. It differs from the standard ap- 
proaches in that it takes the Markov assumption not as an 
approximation in the course of the calculation, but as a 
premise, implemented before tracing out the environment. 
It will be applicable whenever the interaction with the en- 
vironment can reasonably be described in terms of indi- 
vidual interaction events or "collisions" , that is, if one can 
take the environment as consisting of independent (quasi)- 
particles which probe the system each at a time, in the 
sense that both the rate and the effect of an individual 
collision are separately physically meaningful and can be 
formulated microscopically. One may then implement the 
Markov requirement right from the outset by disregarding 
the change of the environmental state after each collision. 
This will be justified if the environment is sufficiently large 
and stationary, and in particular if many different envi- 
ronmental (quasi)-particles are involved so that each has 
much time to carry away and disperse its correlation with 
the system. 

It is clear that the apparatus of time-dependent scat- 
tering theory [7] is predestined for this type of descrip- 
tion. Its microscopically defined S-matrix maps from the 
incoming to the exact outgoing asymptotes of the system- 
environment state without a temporal evolution, and a 
partial trace over the scattered environment yields the sys- 
tem state after a single collision. One would like to write 
the temporal change of the system as the rate of collisions 
multiplied by this change due to an individual scatter- 
ing. The great difficulty with this is that in general also 
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the rate depends on the system state so that a naive im- 
plementation would yield a nonlinear equation. Below, I 
describe how this is circumvented by applying the concept 
of generalized and continuous measurements. The use and 
strength of the method is then demonstrated by deriving 
the master equation for the internal quantum dynamics of 
an immobile system affected by a gaseous environment in 
terms of the multichannel scattering amplitudes. 

Monitoring approach. — My first aim is to argue 
that the system p evolves as d t p = (ih)^ 1 [H, p] + Cp with 
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T + Tt,rV2 [p^p^Y 1 ' 2 
+ Tr cnv (Tr 1 / 2 [p® Pcm ]T 1/2 V) 
— ^ Tr onv (r 1 / 2 T t TT 1 / 2 [p® Pcm ]) 



(i) 



Here H is the Hamiltonian of the isolated system and p cnv 
the reduced single-particle state of the environment. The 
operator T is the nontrivial part of the two-particle S- 
matrix S — I + iT describing the effect of a single collision 
between environmental particle and system. The rate of 
collisions is described by T, a positive operator in the to- 
tal Hilbert space, which determines the probability of a 
collision to occur in a small time interval At, 



Prob(C At |/») 



Attr (r [pOPenv])- 



(2) 



Like the S-matrix, the operator T can in principle be char- 
acterized operationally in independent experiments. Its 
microscopic formulation will in general involve a total scat- 
tering cross section and the current density operator of the 
relative motion (see below). 

To motivate the time evolution (1) we picture the envi- 
ronment as monitoring the system continuously by send- 
ing probe particles which scatter off the system at ran- 
dom times. The state dependent collision rate can now be 
incorporated into the dynamical description by assuming 
that the system is encased by a hypothetical, minimally 
invasive detector with time resolution At. It tells at any 
instant whether a probe particle has passed by and is going 
to scatter off the system. 

The important point to note is that the information 
that a collision will take place changes our description of 
the impinging two-particle state. According to the the- 
ory of generalized measurements [5, 8, 9] the new state 
is the normalized image of a norm-decreasing completely 
positive map A4(-\C At ) in the total Hilbert space satisfy- 
ing ti(M(g\C At )) = Atti(Tg). For an efficient [10] and 
minimally-invasive detector it has the form 



M(g\C At ) 



AtY^gY 1 ' 2 . 



(3) 



The significance of this measurement transformation is to 
imprint our improved knowledge about the incoming two- 
particle wave packet, and it may be viewed as enhancing 



those parts which head towards a collision. In principle, an 
efficient measurement (which introduces no classical noise 
by mapping pure states to pure states) may be given by 
a more general operator, .M(f?|CAt) = M A t£?M^ t as long 
as it satisfies M A( M A i = AtY. The above 'minimally- 
invasive' choice of Ma* is reasonable because a possi- 
ble unitary part Ua* in its general polar decomposition 
Ma* = UAtr^VAt would describe a reversible "back ac- 
tion" which has no physical justification in our case of 
a thought measurement invoked only to account for the 
state dependence of collision probabilities. 

Also the absence of a detection event during At 
changes the state. The corresponding complementary map 
M(-\C A t) satisfies tr (A4(g\C A t)) = l-Attr(Yg) and the 
Kraus representation with time-invariant operators reads 
M(g\C A t) = g-Y 1 / 2 gY 1 ' 2 At. 

We can now form the unconditioned system-probe state 
after time At by allowing for the fact that the detection 
outcomes are not really available. Thus, the infinitcsimally 
evolved state is given by the mixture of the colliding state 
transformed by the S-matrix and the untransformed non- 
colliding one, weighted with the respective probabilities, 



d (At) 



Prob (C At |p)S 



M(g\C At ) 
tr (M(g\C At )) 
M(g\C At ) 



St 



+ Prob ( CA^) tr(M ^ Ai)) 

= sr^v^stAt + e - rVV^At. 

Using the unitarity of S, which implies i (T — T^J 
— T^T, the differential quotient can be written as 

At u 2 y 

- 1 r l/2 er l/2 T t T + i J Rc (T) ) 1-1/201-1/2 



One arrives at (1) by tracing out the environment with 
g = p ® Penv, taking the limit of continuous monitoring 
At — > 0, and adding the generator of the free system evo- 
lution. Thus, the collision rate with its state dependence 
is incorporated by the operators Y 1 / 2 and they may be 
thought of, in a stochastic unravelling of the master equa- 
tion [2, 11-15], as serving to weight each trajectory with 
the rate before it scatters. The operators T describe the 
individual microscopic interaction process without approx- 
imation. Note also that (1) generates a dynamical semi- 
group by construction since M(-\C At ) and M(-\C At ) are 
completely positive. 

To judge whether the trace in (1) yields a useful master 
equation one has to specify system and environment. A 
first application of this general equation can already be 
found in the recent Ref. [16], where it is used to describe 
the motion of a distinguished, freely moving point-particle 
in the presence of a gas. The above discussion thus serves 
to complete the derivation in [16], where a quantum ver- 
sion of the linear Boltzmann equation was obtained which 
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displays all expected limiting properties. In the follow- 
ing, I will demonstrate the use and generality of eq. (1) 
by posing a complementary question, namely, how the the 
internal dynamics of an immobile system gets affected by 
an environment of structureless gas particles. 

Application to an immobile system. If the mo- 
tional system degrees of freedom are disregarded, a single- 
particle S-matrix can be used to describe the (in gen- 
eral inelastic) interaction with the environmental parti- 
cles. The resulting master equation should describe non- 
perturbatively both the coherent and the incoherent pro- 
cesses induced by this coupling. An example would be the 
collisional decay of molecular eigenstates into chiral config- 
urations, or the phonon-induced decoherence of a quantum 
dot. For concreteness, the environment is assumed to be 
an ideal Maxwell gas of density n gas , atomic mass m, and 
single particle state p cnv = (A^ h /r2) exp (— /3p 2 /2m) with 
p the momentum operator, Ath = hyj2-K(3/m the thermal 
wave length, and f2 the normalization volume. 

In the language of scattering theory the free energy 
eigenstates of the non-motional degrees of freedom are 
called channels. In our case, they form a discrete basis 
of the system Hilbert space, and \a) will be used to in- 
dicate internal (and possibly rotational), non-degenerate 
system eigenstates of energy E a . In this channel basis, 
pa/3 = (a\p\(3), the equation of motion (1) takes on the 
form of a general master equation of Lindblad type, 

dtp a p = -p a( i + 2^ PooPo M a p' (4) 

\ E <w E M T \ E p°fh E <7 

c*o 7 fa 7 

with energy shifts e a discussed below and rate coefficients 
Kf° =H Tr c „v (Vr 1 / 2 [|a )(/3o| ® /w] 

(5) 

To calculate these complex quantities we need to spec- 
ify the rate operator Y. In the present case it is naturally 
given in terms of the current density operator j = n gabS p/m 
of the impinging gas particles and the channel-specific to- 
tal scattering cross sections a (p, a) , 

T = V|«)(a|®n gas M CT(P;a ). (6 ) 

a 

Defining the channel operator c = J2 a a \ a )( a \ one can 
thus write r = |j|<r (p, c). 

In principle, T must also involve a projection to the sub- 
space of incoming wave packets, attributing zero collision 
probability to any wave packet located far off the scatter- 
ing center and travelling away from it. This is important 
because such an outgoing state will not remain invariant 
under S. [It may be strongly transformed since the def- 
inition of S involves a backward evolution.] In practice, 



the microscopic definition of T is easier if one takes care 
of the projection separately. This is easily done if p cnv ad- 
mits a convex decomposition into incoming and outgoing 
states. Alternatively, one may dispense with the projec- 
tion by modifying the definition of S so that outgoing wave 
packets are kept invariant (see below). 

Let us now evaluate the rate coefficients M"p by 
using a decomposition of p cnv that permits to separate 
in- and out-wave packets. As shown in [17] the ther- 
mal gas state can be written as a phase space integration 
over projectors onto minimum uncertainty gaussian states 

IVv oPo ) = ^IC ex P (~0 (P _ Po) 2 / 4m ) ko) whose spatial 

extension A t h = fr\j2ir(3/m is determined by an inverse 
temperature > (3, 

Pcnv = y'dpo£(Po)^^IVVoP >WVoP l- (7) 

Here p(p ) — (27rm//3)~ 3//2 cxp(— f3pl/2m) is the 
Maxwell-Boltzmann distribution corresponding to the 
temperature = — /3 _1 , so that by setting a (3 one 
splits up the gas temperature {3~ l into a part determining 
the localization of the |Vv p ) and a part characterizing 
their motion. We choose (3 large and take eventually the 
limit (3 —* oo, (3 — > (3 of very extended wave packets so 
that p approaches the original Maxwell-Boltzmann distri- 
bution p. Inserting (7) into (5) yields 

<r = Jdp Q p(p )J^m a J (ro,Po)- (8) 
Here the phase space function 

<°f°(ro,P ) := |dp<a|<p|TrV 2 |ao>|Vv p ) 

x<M^op |r 1/2 Tt|/?)|p) (9) 

gives the contribution of different phase space regions to 
the rate coefficient M'^ . This permits now to restrict 

the calculation to incoming wave packets. Since the m^f 
are averaged over all available positions in (8) it is natural 
to confine this spatial average at fixed p to a cylinder 
pointing in the direction of p , whose longitudinal support 
A p vanishes at outgoing positions and whose transverse 
base area is given by an average cross section E p . In 
terms of the longitudinal and transverse positions r\\ Po := 
(r ■ Po) Po/Po and r ^P = r - r \\p we have 

J JApo A Po ->£p c ^Po 

x<f°(r\\ Po +r± Po ,P ). (10) 

In order to evaluate Tn"°£° , insert momentum resolutions 
of unity between the T and T operators in (9) and use the 
representation [7] 

(a f \(p f \T\ ai }\ Pi ) = fafa ;}^ Pi) " S(E Pfaf -E Piai ) 

(11) 
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in terms of the mult i- channel scattering amplitude and 
total energies E pa — p 2 /2m + E a . By transforming the 
new integration variables to mid-points and chords one 
obtains a Gaussian function which approaches, for large 
j3, a 5-function in the midpoints. Integrating out the lat- 
ter one finds that the combination of the (^-functions from 
(11) confine the chord integration to a plane perpendicu- 
lar to p . Integrating out the parallel component leads 

to the factor exp (—(3m (E a - E ao — Ep + Ep ) 2 /8poj 

which, again for large (3, can be replaced by 



Aq/3 



1 if E a — E ao — Ep — Ep 
otherwise. 



The resulting expression is independent of T"|| p , 



m a a f°(r ,P ) =X^°5f I dp 



X f aao {p,Po) fho (P'Po) 



\ 



1 + 



a (p£,a ) <J (p ,Po) 



with Pq := p ± p ipo /2. The r|| Po -integration in (10) 
yields an approximate two-dimensional ^-function in p± Po 
so that we obtain 



M 



aoPo _ 
a/3 



~ X a p 



c*o/3o n gas 



TfL 



dp dp p(p )f aao {p,p Q ) 



2 2 

P -Po 

2m 



E„ — E, 



(12) 



provided we identify the average cross section of (10) with 
the geometric mean of the total cross sections of the in- 
volved channels, i.e., S Po = \f a (p ; ceo) o~ (po! A))- More- 
over, the final limit [3 — > oo replaced p by p in (12). 

With the same method one shows that the first term in 
(1) merely modifies the unitary evolution. Its effect is to 
shift the system energies from E a to E a + e a by a thermal 
average of the "forward scattering amplitudes" , 



■ a = -27rft 2 ^£ / dp oM (p ) Re [f aa (pq.Po)] • 

Tfl J 



(13) 



It is reassuring that the explicit expressions (12) and (13) 
can be shown to be equivalent to the more abstract master 
equation by Diimcke [18], obtained in a "low-density limit" 
scaling approach [1,6, 19] for the special case of a factor- 
izing interaction potential, V to t = A® B cnv , and for times 
large compared to all system time scales. The present ap- 
proach thus generalizes this result to arbitrary interaction 



potentials (satisfying asymptotic completeness) and to ar- 
bitrary times as long as they are greater than the duration 
of a single collision. 

It is worth noting that the M*p can as well be ob- 
tained in a more direct, while less solid way if the diagonal 
momentum representation of p CDV is used instead of (7). 
A projection to the incoming wave packets is then hard to 
implement and, as discussed above, the application of S to 
improper momentum states leads to the unwanted trans- 
formation also of its "outgoing components" . As a conse- 
quence, the resulting expression for M°f° is ill-defined, 
involving the square of the 8- functions in (11) and the 
normalization volume 0. This can be healed by noting 
that any consistent modification of S which keeps outgo- 
ing wave packets invariant must conserve the probability 
current. This condition provides a simple rule how to form 
a well-defined expression [17,20], whose multichannel ver- 
sion yields the result (12) immediately for any momentum 
diagonal p cnv . 

The expression for the rate coefficients can be rewrit- 
ten, for isotropic p, in terms of an average over the 
velocity distribution v (v) = Airm 3 v 2 p (mv) and angu- 
lar integrations, which bring about the velocity v out = 
\Jv 2 — 2 (E a — E ao ) jm of the gas particle after a possi- 
bly inelastic collision. For rotationally invariant scattering 
amplitudes, f aao (cos (p,p ) ; E = Po/2m), we have 

/•OO pi 

Kf° =xlf° J Q dv v ( y ) Wout27r J i d (cos 6) 

x f aao (cos 6- - v 2 ) f*pp o (cos 6- -v 2 ) . (14) 

This shows that limiting cases of (4) display the ex- 
pected dynamics. For the populations p aa it reduces to 
a rate equation where the total cross sections a aao (fw 2 ) 
for scattering from channel ceo to a determine the tran- 
sition rates, M^ a a ° = J dvv (v) n gabS v out a aao . In the 
case of purely elastic scattering, on the other hand, i.e., 



M 



aoPo 
a/3 



MapSaagSppQ, the coherences decay exponen- 



'q/3 



tially, d t \p a p\ — — 7 a |f stlc I Pa/3 1 > with a rate determined by 
a difference of scattering amplitudes, 



(15) 



7a| =7F / dv V ( W ) "gas^out / d (COS 6) 



fa 



^v 2 )-fpp( C osO-™v 2 ) 



It shows clearly that the more coherence is lost, in this 
case, the better the scattering environment can distinguish 
between system states \ce) and \(3). 

Conclusions. — In conclusion, a general method of in- 
corporating formal scattering theory into the dynamic de- 
scription of open quantum systems was presented. Based 
on the theory of generalized measurements, it yields com- 
pletely positive master equations which account for the 
environmental interaction in a non-pcrturbative fashion. 
When applied to an immobile system in the presence of 



a gas, it provides a detailed and realistic account of the 
interplay between coherent system dynamics and the (pos- 
sibly much faster) incoherent effects of the environment. 

* * * 
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